The mass-flow error in the Numerical Renormalization Group method 
and the critical behavior of the sub-ohmic spin-boson model 
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We discuss a particular source of error in the Numerical Renormalization Group (NRG) method 
for quantum impurity problems, which is related to a renormalization of impurity parameters due to 
the bath propagator. At any step of the NRG calculation, this renormalization is only partially taken 
into account, leading to systematic variation of the impurity parameters along the flow. This effect 
can cause qualitatively incorrect results when studying quantum critical phenomena, as it leads to 
an implicit variation of the phase transition's control parameter as function of the temperature and 
thus to an unphysical temperature dependence of the order-parameter mass. We demonstrate the 
mass-flow effect for bosonic impurity models with a power law bath spectrum, J(u>) oc uj s , namely 
the dissipative harmonic oscillator and the spin-boson model. We propose an extension of the NRG 
to correct the mass-flow error. Using this, we find unambiguous signatures of a Gaussian critical 
fixed point in the spin-boson model for s < 1/2, consistent with mean-field behavior as expected 
from quantum-to-classical mapping. 



I. INTRODUCTION 



The Numerical Renormalization Group method, 1,2 
originally developed by Wilson 1 for the Kondo model, 
is by now an established technique for the solution of 
general quantum impurity problems. It has been ap- 
plied, e.g., to magnetic atoms in metals, to quantum 
dots and magnetic molecules, and as an impurity solver 
within dynamical mean-field theory. Its generalization 3,4 
to bosonic baths has enabled the treatment of dissipa- 
tive impurity models and those with both bosonic and 
fermionic baths. 5 Quite often, impurity quantum phase 
transitions 6 are in the focus of interest. The strengths 
of NRG in treating such critical phenomena lie in its 
ability to treat arbitrarily small energy scales and in its 
renormalization-group character which allows e.g. for the 
analysis of flow diagrams. 

Recently, conflicting results have been reported about 
the critical behavior of certain impurity models with a 
bosonic bath, in particular the spin-boson and the Ising- 
symmetric Bose-Fermi Kondo model. 7 For a bosonic bath 
with power-law spectral density J(u>) oc lo s , these mod- 
els display a quantum phase transition for < s < 1. 
Statistical-mechanics arguments suggest that this transi- 
tion is in the same universality class as the thermal phase 
transition of the one-dimensional (Id) Ising model with 
l/r 1+s long-range interactions. At issue is the validity 
of this quantum-to-classical correspondence for s < 1/2 
where the Ising model is above its upper-critical dimen- 
sion and displays mean-field behavior. 8,9 Initially two 
of us claimed non-classical behavior with hyperscaling 
in the spin-boson model for s < 1/2, based primar- 
ily on NRG results. 10 These results have been verified 
by others, 11 and extended to the Ising-symmetric Bosc- 
Fermi Kondo model. 5 In contrast, subsequent Quantum 
Monte Carlo (QMC) 12,13 and exact-diagonalization 14 



studies concluded that the critical behavior of the spin- 
boson model for s < 1/2 is classical and of mean-field 
type. We have recently retracted the claim 10 of non- 
classical behavior, because we have realized two different 
sources of error of the NRG which spoil the determination 
of critical exponents. 15 However, other authors continue 
to rely on NRG results in this context. 16,1 

In this paper, we investigate one of the error sources 
of the NRG in more detail, which we have dubbed the 
mass-flow effect. It arises from the NRG algorithm which 
iteratively integrates out the impurity's bath. For a 
particle-hole asymmetric bath, the real part of the bath 
propagator generates a physical shift of impurity param- 
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FIG. 1: Schematic phase diagram of the spin-boson model 
Hsb, as function of temperature T and tunneling strength 
A (keeping the dissipation strength a fixed), showing a path 
with A = A c (thick solid) along which quantum-critical ob- 
servables should be measured. Due to the mass-flow effect, 
a temperature-dependent deviation of the order-parameter 
mass is induced, such that a system located at A c at T — 
follows a path with finite mass, <5A oc T s , at any T > (thick 
dashed). The thin lines represent trajectories with (dashed) 
and without (solid) mass- flow effect for A ^ A c . 
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eters; for models with single-particle tunneling between 
impurity and bath it is simply the energy of the impu- 
rity level which is shifted due to the real part of the 
hybridization function. As a NRG calculation ignores 
the part of the bath spectrum below the current NRG 
scale, there is, at any NRG step, a missing parameter 
shift which is is set by the current NRG scale. Near a 
quantum phase transition, this implies an artificial scale- 
dependent shift of the order-parameter mass. For a NRG 
calculation with model parameter values corresponding 
to the critical point, the system is therefore not located 
at the critical coupling for any finite T, but effectively 
follows a trajectory in the phase diagram as sketched in 
Fig. 1. This spoils the measurement of critical properties 
extracted in a NRG run as function of T. 

Other sources of error within the NRG method are the 
discretization of the continuous bath density of states, 
the truncation of the eigenvalue spectrum in each NRG 
step to the lowest N 3 states, and the truncation of the 
bath Hilbert space in the case of a bosonic bath, where 
only Nb states are taken into account. While the effects of 
discretization and spectrum truncation are well studied 
and understood within the fermionic NRG, 1,2 Hilbcrt- 
space truncation is more serious. Ref. 4 pointed out that 
it precludes a correct representation of the ordered phase 
of the spin-boson model for s < 1 at low energies or tem- 
peratures. Later, it was realized 15 that it also leads to 
incorrect results for the order-parameter exponents (3 and 
5 of the phase transition above the upper-critical dimen- 
sion. In this paper, our focus will be on the mass-flow 
effect; the other errors will be discussed when appropri- 
ate. 

First, we shall demonstrate the mass-flow effect within 
a model of non-interacting bosons, namely the dissipative 
harmonic oscillator, for which all statements can be made 
exact. In this model, the critical point translates into the 
instability point where the renormalizcd impurity energy 
is zero. The mass flow will be shown to lead to qualita- 
tively incorrect results; this problem carries over to inter- 
acting models, like the anharmonic oscillator or the spin- 
boson model, if the critical point is Gaussian (i.e. above 
its upper-critical dimension). Second, we propose an ex- 
tension of the iterative diagonalization scheme to cure 
the mass-flow error. This extension solves the problem 
for the full parameter range of the non-interacting har- 
monic oscillator, while working asymptotically for mod- 
els of interacting bosons. Third, we apply the extended 
NRG algorithm to the spin-boson model. For s < 1/2, 
we find results qualitatively different from those 3 ' 10 of the 
standard NRG implementation: our new results signify 
a flow towards a Gaussian critical fixed point. While the 
truncation of the bosonic Hilbert space precludes calcu- 
lations very close to this Gaussian fixed point, we can 
identify a mean-field power law in the impurity suscepti- 
bility. Taken together, this shows that - as other meth- 
ods - also the NRG predicts that the spin-boson model 
exhibits mean-field behavior for s < 1/2. 

It is worth noting that an observation reminiscent of 



the mass-flow effect has been made in Ref. 16: the criti- 
cal behavior of the classical long-range Ising model with 
s < 1/2 was found to change from mean-ficld-like to 
hyperscaling-like upon artificially truncating the "wind- 
ing" of the long-range interaction (i.e. upon violating the 
periodic boundary conditions in imaginary time). This 
finding underscores that mean-field critical behavior in 
the long-range models under consideration can be easily 
spoiled by algorithmic errors. Note, however, that we 
disagree with the interpretation regarding the quantum- 
to-classical correspondence given in Ref. 16, see below. 

A. Outline 

The remainder of the paper is organized as follows: In 
Sec. II the model Hamiltonians are introduced. Sec. Ill 
explains how the mass-flow effect arises from the iterative 
diagonalization of the Wilson chain. The dissipative har- 
monic oscillator is subject of Sec. IV, where the mass-flow 
error in the susceptibility is demonstrated analytically. 
This knowledge is used in Sec. V to propose a modifica- 
tion of the iterative-diagonalization scheme, designed to 
cure the mass-flow error. Finally, the modified NRG algo- 
rithm is applied to the spin-boson model in Sec. VI. The 
NRG flow is discussed separately for s > 1/2 and s < 1/2 
and compared to the results from standard NRG. The re- 
sults are interpreted in terms of a Gaussian critical fixed 
point for s < 1/2. Conclusions close the paper. Various 
details, including a discussion of the mass-flow effect in 
fermionic impurity models with particle-hole asymmetry, 
are relegated to the appendices. 

II. MODELS 

The mass-flow effect can be most easily demonstrated 
using impurity models of non-interacting particles. We 
shall consider the dissipative harmonic oscillator, with 
the Hamiltonian 

Hdho = £!a t a+ |(a + « f ) 

+ l^Kia + a^ibi + b^ + ^Uibtbi, (1) 

i i 

where 17 > is the bare "impurity" oscillator frequency, e 
is a field conjugate to the oscillator position, and the Wj > 
are the frequencies of the bath oscillators. The bath is 
completely specified by its propagator at the "impurity" 
location 

I» = V ^ (2) 

with the spectral density 

J(oj) = -lmT(oj) = 71^2^6 (lu-lu^ . (3) 
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Universal properties of impurity phase transitions are 
determined by the behavior of the low-energy part of the 
bath spectrum J(w). Discarding high-energy details, the 
common paramctrization is 



2naujl~ s uj s 



< LU < LU C 



s > -1 (4) 



where the dimensionless parameter a characterizes the 
dissipation strength, and oj c is a cutoff energy. The value 
s = 1 represents the case of ohmic dissipation. 

The dissipative oscillator with a power-law bath spec- 
trum is known to become unstable at large dissipation: 21 
The coupling to the bath renormalizes the oscillator fre- 
quency Q downwards, which becomes zero at some a c . 
Hence, the behavior of the model is not well-defined for 
a > a c . 

The system at large dissipation may be stabilized 
by adding a local repulsive interaction to 7Ydho- A 
symmetry-broken phase can emerge, with "condensa- 
tion" of the a bosons. Two possible routes are 



Hdao = Hdho + Un a (n a - 1), 



and 



n' 



DAO 



Hdho + u(a ■ 



n„ = a) a 



^4 



(5) 



(6) 



The latter, 7ip A o, can be understood as a local <f> im- 
purity. On the other hand, 7Ydao m the limit U = oo 
becomes equivalent to the standard spin-boson model 

Wsb = - ^ ^ + + y Yl \{k + b\) + Y,u4bi (7) 



where a z = ±1 are the local impurity states, and fl is 
the tunneling rate. The equivalence is seen by identifying 
the remaining oscillator states |0) and |1) in 7Ydao with 
the states (| |) ± | l))/V%- In all three models (5,6,7), 
the ordered phase at large dissipation breaks an Ising 
symmetry, a <-* —a (or a z <-* — er 2 ), b <-> —6, and is 
associated with a non-zero expectation value (a + a^) (or 

Universality arguments suggest that the critical prop- 
erties of the phase transitions are identical in the three 
models and coincide with those of a classical Ising chain 
with l/r 1+s interactions. This quantum-to-classical cor- 
respondence trivially holds for "H'^^q in Eq. (6), as its 
imaginary-time path integral representation at T = is 
identical to the continuum limit of the one-dimensional 
Ising model (i.e. a scalar 4> 4 theory). 9 For s < 1/2, the 
critical behavior is Gaussian and mean-held like, with 
the quartic interaction being dangerously irrelevant at 
criticality. 

The quantum phase transition in the spin-boson model 
has been extensively studied: While the ohmic case, 
s = 1, has long been known to display a Kosterlitz- 
Thouless transition, 22 the sub-ohmic case has only been 
investigated more recently 3,10 ' 23 For < s < 1, a con- 
tinuous quantum phase transition emerges, with critical 



exponents depending on s. While there is consensus that, 
for 1/2 < s < 1, those exponents are identical to the ones 
of the corresponding Id Ising model, a debate is centered 
around the issue of whether or not this continues to hold 
for < s < 1/2 where the Ising model displays mean- 
held behavior. Alternatively, non- mean-held exponents 
obeying hyperscaling have been proposed on the basis of 
NRG calculations 10,11 , and also carried over to the Ising- 
symmetric Bose-Fermi Kondo model. 5,7,16 In particular, 
NRG has been used to calculate the local susceptibility \ 
at the critical coupling as function of temperature, which 
was found to follow a power law \ cx T~ x with x = s. 
In contrast, mean-held behavior implies 8 x = 1/2, which 
has indeed been found e.g. using QMC simulations. 12 

We shall argue here, expanding on our previous note, 15 
that the proposals of non-classical behavior are erroneous 
for the spin-boson model and questionable for the Ising- 
symmetric Bose-Fermi Kondo model. For the former, we 
show that the critical behavior instead is of mean-held 
type, consistent with numerical studies using QMC and 
exact-diagonalization methods. 12-14 



III. WILSON CHAIN AND MASS FLOW 

Within the NRG algorithm, the bath is represented by 
a semi-infinite ("Wilson") chain, Fig. 2, such that the 
local density of states at the hrst site of this chain is a 
discrete approximation to the bath density of states. 1,2 
Due to the logarithmic discretization, the site energies 
e n and hopping matrix elements t n decay exponentially 
along the chain according to w c A~ n+1 , where A is the 
discretization parameter. 

Let us denote by Tl n the Hamiltonian of impurity plus 
n sites of the Wilson chain, and by T n (uj) the propaga- 
tor at the impurity site of this n-site bath. Then, Tioa is 
the discretized version of the original problem. During 
the NRG run, Hoc is diagonalized iteratively: First, Hi 
is diagonalized and the lowest N s eigenstates are kept. 
Then, the next bath site is added to form TI2, the new 
system is diagonalized, and again the lowest N s eigen- 
states are kept (which are approximations to the lowest 
states of 7^2)- As the characteristic energy scale of the 
low-lying part of the eigenvalue spectrum decreases by 
a factor of A in each step, this process is repeated un- 
til the desired lowest energy is reached. Temperature- 
dependent thermodynamic observables at a temperature 




FIG. 2: Structure of the NRG Hamiltonian, with the bath 
represented by a semi-infinite Wilson chain, with bath op- 
erators b' n . The boxes indicate the iterative diagonalization 
scheme. 
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FIG. 3: a) Zero-frequency real part of the bath propagator 
at the initial site of the Wilson chain, — ReG^(o; = 0), as 
function of the chain length n, for different bath exponents s 
and discretization parameters A, all with uj c = 1. b) The 
piece of the real part which is missing at chain length n, 
— Re (G^o — G°), for the same parameters, plotted as func- 
tion of the characteristic energy scale of the Wilson chain, 



T n — ui c A~ n+1 1 (3 are typically calculated via a thermal 
average taken from the eigenstates at NRG step n. Here, 
[3 is a parameter of order unity which is often chosen as 
(3=1. 

The iterative diagonalization procedure implies that, 
at NRG step n, the chain sites n + 1, n + 2, . . .have 
not yet been taken into account, i.e., the effect of those 
sites does not enter thermodynamic observables at tem- 
perature T n . Typically, this is a reasonable approxima- 
tion, as the spectral density of the missing part of the 
chain, Iin^oo — T n )(u>), has contributions at energies be- 
low LO c hr n only. 

However, the missing chain also implies a missing con- 
tribution to the real part of the bath propagator. This 
can be easily estimated: For a power-law bath spectrum, 
Eq. 4, the zero-frequency real part Re(r oo — T n )(uj = 0) 
is generated by frequencies < u> < uj c A~ n and scales 
as uj c A~ ns , i.e., up to numerical factors it scales as the 
NRG energy scale T n to the power s. As we will show 
below, this missing real part implies a flow of the order- 
parameter mass and can spoil the analysis of critical phe- 
nomena. 

To support the above estimate, we calculate the local 
Green's function G° at the initial site, b' 1: of the Wilson 
chain (which is proportional to T n (u))) for different chain 
lengths n. To this end, we numerically diagonalize the 
single-particle problem corresponding to a Wilson chain 
with parameters e n and t n chosen to represent a power- 
law bath spectrum, Eq. (4), as in the NRG. 24 

Explicit results for ReG°(o; = 0) are shown in Fig. 3a. 
As expected, Re G° approaches a finite (negative) value 
as n — > oo, which depends on both s and A. The missing 
real part Re (G^ — G°) is shown in panel b and scales as 
T*, with a prefactor which depends on s, but only weakly 
on A. 



IV. DISSIPATIVE HARMONIC OSCILLATOR 

We shall discuss how the mass-flow effect influences 
observables for the simplest model, the dissipative har- 
monic oscillator (1). It is important to distinguish the 
various methods to calculate observables in this non- 
interacting model: (i) For a continuous power-law spec- 
trum, a number of quantities can be calculated analyti- 
cally, (ii) For a discretized bath, represented by a semi- 
infinite Wilson chain, the single-particle problem can be 
solved by exact diagonalization for long chains, (iii) As 
in NRG, one may use a truncated Wilson chain with 
temperature-dependent length, and again diagonalize the 
single-particle problem, (iv) A true NRG calculation can 
be performed, which treats the full many-body problem. 
Here, we shall mainly be interested in comparing the re- 
sults of (ii) and (iii), which allows to assess the mass-flow 
error. In contrast, the difference between (i) and (ii) can 
be used to quantify the discretization error, while the 
difference between (iii) and (iv) is due to spectrum and 
Hilbert-space truncation of NRG. 

The most interesting observable is the susceptibility 
associated with the oscillator position, defined according 
to 



X = d(a + a))/de, 



(8) 



which is the analogue of d(a z )/de in the spin-boson 
model. Importantly, \ is given by a single-particle prop- 
agator, x = —G x , with 



G x (uj) = ((a + a f ; a + a f ) 



(9) 



2n 



Q 2 -fi[T(w) + r(-w)]/2' 



note the factors of 1/2 in Eq. (1). This equation shows 
that the dissipative oscillator is unstable at and be- 
yond the "resonance" which occurs at some dissipation 
strength a c , defined by Q + Rer(o; = 0) = 0. For a < a c , 
all eigenenergies of the system are positive, whereas the 
lowest one turns negative for a > a c . Thus, a c corre- 
sponds to a singularity of the dissipative harmonic oscil- 
lator, separating the stable from the unstable regime. 

Returning to the susceptibility x, its static limit eval- 
uates to 



X 



Q + RcT(uj = 0) 



(10) 



which is seen to be temperature-independent and only 
determined by and the real part of the bath propaga- 
tor. Consequently, there is a strong mass-flow effect, as 
the renormalized oscillator frequency in the denominator 
of x reads f2 r = O + Rer„(w = 0) for a n-site chain. 
This is illustrated in Fig. 4 where we show the suscep- 
tibility as function of temperature, calculated using ei- 
ther a long Wilson chain for all T or an n-site Wilson 
chain at temperature T ra , i.e., using methods (ii) and (iii) 
described above. Most importantly, x calculated from 
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T 

FIG. 4: Mass-flow error for the dissipative harmonic oscil- 
lator: the graph shows the susceptibility x> Eq. (8), calcu- 
lated at temperature T n by diagonalizing the single-particle 
problem of a truncated n-site Wilson chain. Parameters are 
Q = 1, s = 0.4, ujc = 1, A = 4 where a c = 0.2284682. The 
different curves are for a = 0.1, 0.22, 0.228, a c (from bottom 
to top). The dashed lines denote the (exact) temperature- 
independent x f° r the a < a c , obtained from the semi-infinite 
Wilson chain. As an aside, we note that the critical coupling 
for a continuum bath spectrum is given by a c — sQ/(2u c ). 
This evaluates to 0.2, showing that the discretization error is 
around 15%. 



the truncated Wilson chain is temperature-dependent, in 
contrast to the exact result. For a < a c , the exact re- 
sult is approached at low T . The error is most drastic at 
resonance, a = a c . There, Xoxact = oo (i.e. the system 
is unstable), whereas the calculation using a truncated 
Wilson chain gives \ °c T~ s . Physicswise, the mass- flow 
effect introduces a finite and temperature-dependent os- 
cillator frequency Q r oc T s , thus artificially stabilizing 
the system at a c . Naturally, the same result is found 
using a full NRG calculation. 

As the dissipative harmonic oscillator represents the 
fixed-point Hamiltonian of the Gaussian critical point 
of, e.g., the anharmonic oscillator Hqao m -^q. (6), it 
is straightforward to discuss the mass-flow effect there. 
The renormalized f2 r can be identified with the order- 
parameter mass, and x( w = 0) = l/O r . Along the flow 
towards the Gaussian fixed point, the irrelevant interac- 
tion u leads to an order-parameter mass oc T 1 / 2 for s < 
1/2, and the physical susceptibility follows x ^ T -1 / 2 
(Ref. 8). However, the artificial mass oc T s caused by 
the mass-flow effect dominates the physical mass at low 
T, leading again to the unphysical result \ ^ T~~ s . As 
this coincides with the physical result for an interact- 
ing critical fixed point with hyperscaling, the unphysical 
result from the mass-flow effect could be mistaken as a 
signature of interacting quantum criticality. 

We should emphasize that a renormalization-group 
scheme which successively integrates out the impurity's 
bath is perfectly valid. However, it requires that the cal- 
culation of observables at some scale T accounts for the 
remaining part of the bath. The latter is not the case in 
the iterative diagonalization scheme of standard NRG. 



V. CURE OF MASS FLOW 

The mass-flow error arises from the missing real part 
of the bath propagator, HeT n (uj = 0), which, for every 
step of the iterative diagonalization, is simply a num- 
ber. Ideally, a general algorithmic solution of the mass- 
flow problem would directly correct Rer„. However, this 
is limited by Kramer-Kronig relations, and we have not 
found a manageable implementation of this idea. 

In the following, we shall instead make use of physics 
arguments in order to (approximately) correct the mass- 
flow error. For the harmonic oscillator, Rer„ directly 
renormalizes the oscillator's energy, while things are con- 
ceptually more complicated for interacting models (like 
the spin-boson or Bose- Fermi Kondo models). There- 
fore, we shall separately discuss the non-interacting and 
interacting cases in the following. 

A. Dissipative harmonic oscillator 

A simple recipe can be used to correct the mass-flow er- 
ror when diagonalizing a finite-length chain correspond- 
ing to the harmonic-oscillator Ti n . We define a Hamilto- 
nian piece JC n by 

K n = TZRe{T 00 -T n )(u; = 0) (11) 

with 1Z = a^a. As a result, H n + IC n has the correct mass 
term, i.e., the correct renormalized oscillator frequency, 
for any n, and diagonalizing H n + K n instead of H n in 
step n removes the mass-flow problem. One obtains the 
correct result for \ : thanks to /C„, the denominator of 
X in Eq. (10) is replaced by ft + Re(r oo - T n ){w = 0) + 
Rer„ (uj = 0) which is the exact result for the semi- infinite 
Wilson chain. 



B. NRG implementation 

A mass-flow correction via K, n can be straightforwardly 
implemented into the iterative diagonalization scheme of 
the NRG method. The modified NRG algorithm (dubbed 
NRG* in the following) works as follows: (i) Initially, one 
diagonalizes H" — TL\ + }C\. In addition to the usual 
observables, the matrix elements of the operator 1Z are 
stored as well. Then, the following steps are repeated: 
(ii) From the lowest N s states of the solution of NRG step 
n and the states of the impurity site n+1, one constructs 
Ti' n+1 . In contrast to TL n +\, the operator Ti.' n+ i contains a 
mass-flow correction from the previous steps, (iii) Using 
the matrix elements of 1Z, one constructs = Ti.' n , 1 + 

K n +i — IC n . (iv) One diagonalizes 7i" +1 and re-calculates 
the matrix elements of the desired observables and of 1Z. 

The correction of the mass-flow error, contained in 
steps (i) and (iii) which differ from the usual NRG al- 
gorithm, is implemented such that the frequency shifts 
cancel in the limit n — > oo. Hence, runs of NRG and 
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NRG* with the same model parameters should target 
the same point in the phase diagram as T — > (although 
their finite-temperature trajectories are different, Fig. 1). 
However, this is only true in the absence of spectrum 
truncation. For finite N s , the cancellation is only ap- 
proximate, i.e., there will be a small (but unimportant) 
parameter shift due to the mass-flow correction. 

C. Beyond non-interacting bosons 

Being interested in extracting critical properties, we 
identify the mass-flow effect as a scale-dependent shift of 
the order-parameter mass. This suggests that the mass 
flow can be corrected by an appropriate shift in the phase 
transition's control parameter - this is simply a gener- 
alization of Eq. (11) where K n shifts the oscillator fre- 
quency. We thus propose to employ a correction of the 
form 

]c n = nn Rc(r oc - r n )(w = o) (12) 

where 1Z is now a (local) operator which can be used 
to tune the phase transition, e.g., the tunneling term a x 
in the spin-boson model or the Kondo coupling term in 
a Bose-Fermi Kondo model. Importantly, the required 
shift will no longer be identical to Re(Foo — T n ). This is 
already clear for the dissipative anharmonic oscillators, 
Eqs. (5) and (6), where the quartic interaction will renor- 
malize both the oscillator frequency and also its shift due 
to the bath, but in a different fashion. Hence, we have 
introduced the non-universal prefactor n which we intend 
to determine by physical criteria. 

Two issues require special consideration: (a) Is the 
linear relation between the required shift in the control 
parameter and the missing real part of T, which is implied 
by Eq. (12), justified? (b) How can one determine the 
prefactor k? 

The simplest argument for (a) is as follows: The 
phase transition's control parameter (equivalently, the 
distance to criticality or the bare order-parameter mass) 
depends on both the prefactor of 1Z and the real part of 
Rer(w = 0). Both dependencies have a regular Taylor ex- 
pansion at a given point in parameter space, hence, the 
leading terms are linear. As Rer changes by a known 
amount in every step of the iterative diagonalization due 
to the mass-flow effect, this can be compensated by a 
change in the prefactor of 1Z which proportional to this 
amount, i.e., a change of the form K n+ \ — JC n with some 
fixed k. This argument only relies on the Taylor expan- 
sion and is thus asymptotically correct for small changes 
in Rer, i.e., for T — > 0. (For a given model, like the 
anharmonic oscillator (6), one can check the linear be- 
havior by an explicit perturbative calculation.) Physi- 
cally, it is clear that the linear term of the expansion will 
capture the correct behavior in the vicinity of a given 
renormalization-group fixed point, i.e., the required n de- 
pends on the fixed point of interest (and on non-universal 
high-energy details). Note, however, that the procedure 



is more general than these considerations suggest: As 
both the Gaussian critical fixed point and the delocalized 
fixed point are asymptotically non-interacting, a fixed k 
can be used to capture the entire crossover from the quan- 
tum critical to the delocalized regime in this case. 

Question (b) will be discussed for different types of 
critical fixed points in turn. We shall use the language 
of the dissipative anharmonic oscillator (6), where the 
critical theory is known. 8,9 

1. Gaussian critical fixed points 

A Gaussian critical fixed point, realized for s < 1/2, 
provides a simple criterion to find the correct value kq 
of the correction parameter k, namely the temperature 
dependence of the order parameter mass. As emphasized 
in Sec. IV, the artificial mass generated by the mass-flow 
effect follows T s , while the physical mass scales as T 1 / 2 . 
Thus, in general the mass at the critical coupling will be 
given by Ai(k — k)T s + A2T 1 / 2 where Ai,2 are prefac- 
tors. For k < n (undercompensation), the positive T s 
term will always dominate at low T and mimic hyper- 
scaling properties. For k > kq (overcompensation), the 
mass will become negative at low T, i.e., the flow will be 
towards the localized phase. An intermediate flow inside 
the localized phase will even occur if couplings are cho- 
sen to be slightly in the delocalized phase: for A > A c 
or a < a c the system flows from critical to localized and 
then back to delocalized upon lowering T, accompanied 
by a non-monotonic behavior of \. 

This suggests the following simple recipe to determine 
kq: Start with large k such that non-monotonic flows are 
seen near the critical coupling. Decrease k until those 
disappear and the susceptibility follows a power law dif- 
ferent from hyperscaling at the critical coupling down to 
the lowest accessible temperatures. If k is decreased too 
far, then \ T~ s is recovered. Hence, a clear signature 
of Gaussian criticality is a qualitatively different behav- 
ior in x f° r small and large k. In Sec. VI and App. A, 
we shall demonstrate this for the spin-boson model at 
s < 1/2. 

2. Interacting critical fixed points 

In the case of an interacting critical fixed point, re- 
alized for s > 1/2, hyperscaling is fulfilled on physical 
grounds. Hence, the mass will invariably scale as T s at 
the critical coupling, both for k < k and k > k - This 
simply reflects the fact that the mass-flow effect does not 
introduce qualitative (but only quantitative) errors here, 
in contrast to the case of Gaussian criticality. Hence, the 
behavior in the quantum-critical regime does not provide 
a sharp physical criterion to determine kq. We conclude 
that a clear signature of true interacting criticality is an 
insensitivity to the value of k of the qualitative critical 
behavior. 



For the spin-boson model at s = 1, a comparison of ob- 
servables to those from other solutions like Bethe Ansatz 
or bosonization could be used to determine kq (for ei- 
ther the localized or the delocalized phase). As s = 1 
plays the role of a lower-critical dimension, we have not 
followed this route further. 



VI. SPIN-BOSON MODEL 

We now apply the modified NRG* algorithm, which 
includes the mass-flow correction (12), to the spin-boson 
model. Note that we will make no a-priori assumptions 
on the nature of the critical fixed points, but instead 
apply the strategies outlined in Sec. VC to determine 
the optimal kq within the NRG* algorithm. 
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A. Flow diagrams 

We have studied the flow diagrams for various values of 
the bath exponent < s < 1 and the mass-flow correction 
parameter k. While a detailed set of data is displayed in 
App. A, the main conclusion is that for s < 1/2 the 
flow changes qualitatively as k is varied, while this is not 
the case for s > 1/2. The former fact can be used to 
determine kq for s < 1/2, while a rough estimate of kq 
for s > 1/2 may be obtained from an extrapolation of 
k (s). 

Doing so, we obtain the flow diagrams from the mass- 
flow corrected NRG* algorithm, which represent a central 
result of this paper. Those are shown in Figs. 5 and 6 
for s — 0.4 and 0.6, respectively, together with the flow 
diagrams from standard NRG. The latter are similar to 
the ones shown in earlier papers. 3 ' 4 

Let us start the discussion with Fig. 5a, displaying the 
standard NRG flow for s = 0.4 near the critical coupling 
strength. For a < a c (left) the flow reaches the delocal- 
ized fixed point, whereas it is directed towards the local- 
ized fixed point for a > a c (right); note that the latter is 
not correctly described due to Hilbert-space truncation. 4 
The flow at a c (dashed) shows a different NRG fixed 
point, which has been identified with the critical fixed 
point. For both a < a c and a > a c this level structure 
is visible at intermediate stages of the flow, before the 
system departs towards one of the stable fixed points - 
this crossover is usually identified with the quantum crit- 
ical crossover scale T* above which the system is critical. 
Now consider the flow of NRG*, Fig. 5b, which includes 
the mass-flow correction of Sec. V. While the asymptotic 
fixed points for both a < a c and a > a c are identical, the 
flow near criticality is strikingly different. In particular, 
no stable level pattern emerges, possibly corresponding 
to a critical NRG fixed point. Instead, all levels appear 
to converge toward zero energy before the critical regime 
is left. Note that the critical flow cannot be followed to 
large n (the system is always localized or delocalized for 
n>20). 



FIG. 5: NRG flow showing the 10 lowest levels for the spin- 
boson model at s — 0.4 obtained from a) standard NRG 
without mass-flow correction and b) NRG* with k, — 0.5 
mass-flow correction. The NRG parameters are A = 4, 
Nb = 12, N s — 40, while the remaining model parameters are 
A = u c = 1. In a), a = 0.3555 (left) and 0.3557 (right), while 
in b) ol = 0.35739 (left) and a = 0.35745 (right). The dashed 
lines show the (near-) critical flow, in a) a c = 0.3555842 and 
in b) a c ~ 0.357992. The flows near criticality in a) and b) 
are qualitatively different. No critical fixed-point structure 
emerges in b), instead the level spacing decreases along the 
flow - this is a signature of the flow towards a Gaussian fixed 
point. (The critical flow in b) cannot be followed beyond 
n ~ 20 due to Hilbert-space truncation errors.) 




5 10 15 20 5 10 15 20 

n n 

FIG. 6: As in Fig. 5, but for s = 0.6 and k = 0.8 in b). 
In a), a = 0.6702 (left), 0.6704 (right), a c = 0.6702965 
(dashed), while in b) a = 0.6723 (left), a = 0.6725 (right), 
and q c = 0.67240901 (dashed). Here, the flows in a) and 
b) are qualitatively similar, with the level structure at the 
critical fixed point in b) slightly deviating from that in a). 
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In Fig. 6, the same comparison of flow diagrams is 
given for s = 0.6. Here, no qualitative difference between 
the flows without and with mass-flow correction is seen. 
A stable level pattern is visible near criticality in both 
cases, but the level energies differ slightly in Fig. 6a and 
b. We found this behavior to be generic for 1/2 < s < 1, 
while the absence of a critical NRG fixed point as in 
Fig. 5b is characteristic for all < s < 1/2, if k is chosen 
according to the criteria in Sec. VC. 



It is straightforward to discuss what would be expected 
for a quantum phase transition above its upper-critical 
dimension. The Gaussian fixed point features free mass- 
less bosons, and interactions are required to stabilize the 
system at T > 0. Those are dangerously irrelevant and 
flow to zero in the critical regime, with a scaling dimen- 
sion which is small near the upper-critical dimension. 
Translated into a many-body spectrum, this implies that 
at the Gaussian critical fixed point the spectrum con- 
sists of an infinite number of degenerate levels at zero 
energy, while the flow towards the critical fixed point 
is characterized by the level spacing flowing to zero as 
n — > co. The latter is precisely what is seen in Fig. 5b. 
It is also clear that within NRG* the fixed point itself 
can never be reached, because with decreasing interac- 
tions (i.e. decreasing level spacing) the error introduced 
by the Hilbert-space truncation becomes more and more 
serious (i.e. bosonic occupation numbers become large). 
This implies that small values of T* cannot be reached 
(as the system always flows to either the localized or de- 
localized phase below some T^ in ) which also limits the 
precision with which we can determine a c . 



A few remarks are in order: (i) During the flow to- 
wards the Gaussian fixed point, Fig. 5b, the rate of de- 
crease in level spacing as function of n depends strongly 
on s, i.e., the level spacing decays faster with smaller s, 
qualitatively consistent with the scaling dimension of the 
interaction u being 8 (2s— 1). Correspondingly, the critical 
flows breaks down earlier for smaller s. (ii) The value of 
the critical coupling a c differs between NRG and NRG* . 
As discussed above, this is a result of spectrum trunca- 
tion within NRG*. We have checked that the difference 
decreases with increasing N s . Further the difference is 
larger for smaller s, which follows from the mass-flow 
error itself being larger for smaller s, see Fig. 3b. 



We conclude that the critical behavior of the spin- 
boson model for s < 1/2 is Gaussian. The stable critical 
fixed point in Fig. 5a is then an artifact of the mass-flow 
error, where the system follows the thick solid trajectory 
in Fig. 1. In contrast, for s > 1/2 the critical theory of 
the spin-boson model is interacting. These conclusions 
are supported by the analysis of xC^Oi see next subsec- 
tion. 
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FIG. 7: Local susceptibility \ °f the spin-boson model at 
s = 0.4 obtained from a) standard NRG without mass-flow 
correction and b) NRG* with k — 0.5 mass-flow correction. 
The a values are a) 0.355, 0.3555, 0.35557, 0.3555842 = a c , 
0.3556, b) 0.357, 0.3573, 0.35739, 0.3573992 w a c , 0.35741, 
with the critical (non-critical) \ plotted with thick (thin) 
lines. The other parameters are as in Fig. 5. The dashed 
(dash-dotted) lines show power laws with T~ s (T~ 1//2 ) as ref- 
erence. The squares in a) show the x of a harmonic oscilla- 
tor, calculated with a truncated Wilson chain corresponding 
to the Qf c of NRG and f2 tuned to resonance. The critical \ 
in panel b) is seen to approach mean-field behavior at low T, 

B. Susceptibility 

We continue with NRG results for the order-parameter 
susceptibility 

X = d(a z )/de (13) 

of the spin-boson model. We will focus on the power-law 
behavior x(T) cx T~ x in the quantum-critical regime. 

For both s = 0.4 and s = 0.3, data from both standard 
NRG and NRG* are shown in Figs. 7 and 8, respectively. 
As reported before, x = s is obtained from NRG, Fig. 7a 
and Fig. 8a, while the correct result near a Gaussian fixed 
point is x — 1/2. It should be noted that this T -1 / 2 
power law requires the renormalized quartic interaction 
to be small. However, once the effective interaction be- 
comes small in the numerics, the NRG* algorithm breaks 
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10" 12 10" 10 10~ 8 10" 6 10" 4 10~ 2 
T 

FIG. 8: Same as Fig. 7, but for s = 0.3 and k = 0.36 in 
b). The a values are a) 0.23, 0.238, 0.2388, 0.238885 = a c , 
0.239, b) 0.24, 0.2408, 0.24083, 0.24085 w a c , 0.2409. Note 
that T^, in « 10 -11 , below which the near-critical curves are 
affected by Hilbert-space truncation. 



down due to Hilbert-space truncation. Thus, the weakly 
interacting Gaussian critical regime cannot be reached, 
and we cannot expect to see an asymptotic T~~ x / 2 sus- 
ceptibility power law. For our parameter values, the 
truncation- induced lower cutoff scale, T* lin , for the crit- 
ical regime is 0(l(r 13 ) for s = 0.4 and ©(lCT 11 ) for 
s = 0.3. Notably, the NRG* results in Figs. 7b and 
Fig. 8b do follow T -1 / 2 over two to three decades in tem- 
perature above T^ in , while T~ s is never seen. 

We are again forced to conclude that the T~ s behavior 
in standard NRG, Figs. 7a and 8a, is an artifact of the 
mass-flow error. To support this, we also show the sus- 
ceptibility of the dissipative harmonic oscillator model, 
Eq. (10), calculated using a truncated Wilson chain with 
the same chain parameters as in the NRG run for the 
spin-boson model and f2 tuned to resonance. As ex- 
plained in Sec. IV, this model has x — °°j but a finite 
X results exclusively from the mass-flow error. Remark- 
ably, this x matches the x from NRG for the spin-boson 
model at low temperatures to an accuracy of better than 
15% - this is consistent with the assertion that the latter 
reflects the physics of a Gaussian fixed point artificially 
stabilized by the mass-flow effect. 




10" 12 10" 10 10" 8 10" 6 10" 4 10" 2 
T 

FIG. 9: Same as Fig. 7, but for s = 0.6 and k = 0.8 in b). 
The a values are a) 0.67, 0.6702, 0.67029, 0.67029648 = a c , 
0.6703, b) 0.67, 0.672, 0.6724, 0.67240901 = a c , 0.67242. The 
dashed line shows a power law with T~ s . 

Susceptibility data for s — 0.6 are shown in Fig. 9. 
Both NRG and NRG* yield a power law with x = 
s, albeit with prefactors differing by 10%. Here, the 
harmonic-oscillator x (with mass flow) and the NRG x 
do not match, but instead differ by roughly a factor of 
1.9. All this is consistent with a true interacting fixed 
point. 



C. Other observables and exponents 

From the discussion, it is obvious that other observ- 
ables at criticality will suffer the mass-flow error simi- 
lar to x(7 1 ). This applies to thermodynamic quantities 
including entropy and specific heat, but also to zero- 
temperature dynamic quantities, like the susceptibility 
x(^). While the latter is defined from the ground state, 
the corresponding NRG evaluation is in fact done during 
the flow. 2 Hence, x( Cl ') f° r s < 1/2 is potentially incor- 
rect as well. However, x(w) can be proven to follow uj~ s 
for all s, irrespective of whether the fixed point is Gaus- 
sian or interacting, 9,25 such that the mass-flow effect only 
introduces quantitative deviations. 

Off-critical properties are to leading order not affected 
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by the mass-flow error, because the artificial mass van- 
ishes as T — > while the physical mass remains finite. 
However, subleading corrections are subject to the mass- 
flow error. 

The NRG calculations of Ref. 10 did not only find the 
critical exponent x to deviate from its mean-field value 
for s < 1/2, but also the order-parameter exponents (3 
and 8. As discussed in detail in Ref. 15, this incorrect 
result is due to a different failure of the bosonic NRG, 
namely the fact that the Hilbert-space truncation pre- 
vents an asymptotically correct representation of the lo- 
calized fixed point for s < 1. For mean-field criticality, 
(3 and 8 are not properties of the critical fixed point, but 
instead of the flow towards the localized fixed point. As 
the latter suffers from the Hilbert-space truncation, (3 
and 8 are unreliable. However, large values of Nb can be 
used to uncover the physical power laws at intermediate 
scales (which are of mean-field type for s < 1/2), before 
truncation effects set in. 15 



VII. OTHER MODELS 

For both versions of the anharmonic oscillator, Eqs. (5) 
and (6), we have obtained results which are qualitatively 
similar to those for the spin-boson model. In particular, 
the standard NRG exhibits signatures of an interacting 
critical fixed point for all s, as in Figs. 5a and 6a. For 
^dao> this result is obviously incorrect for s < 1/2, due 
to its equivalence to a local (f> 4 theory. Accordingly, the 
mass-flow corrected NRG* algorithm yields Gaussian be- 
havior in both models for s < 1/2. Hence, all models 
(5,6,7) belong to the same universality class and follow 
the quantum-to-classical correspondence. 

We have also investigated the mass-flow effect for 
particle-hole asymmetric fermionic impurity models. 
While generically present, its effects on observables turn 
out to be tiny, for details see App. B. 

Finally, a remark on symmetries and the quantum- 
to-classical correspondence is in order: While all cases 
discussed so far feature Ising-symmetric critical degrees 
of freedom, impurity spin models with higher symmetry 
[e.g. SU(2)] have been discussed extensively in the litera- 
ture as well. Here, a direct quantum-to-classical mapping 
(via a re-interpretation of the Trotter-discretized action 
of the quantum model after integrating out the bath) is 
usually not possible, due to the impurity spin's Berry 
phase. Indeed, the so-called Bose-Kondo model with 
SU(2) symmetry exhibits a stable intermediate-coupling 
fixed point 26 (unlike any classical Id spin model), and 
the SU(A)-symmetric Bose-Fermi Kondo model has been 
shown to display a quantum critical point with hyperscal- 
ing for all s. 27 



VIII. CONCLUSIONS 

In this paper, we have investigated a source of error in 
Wilson's NRG method which had received little attention 
before. This mass-flow error is inherent to the iterative 
diagonalization scheme of NRG which neglects the low- 
energy part of the bath when calculating observables. 

We have traced the mass-flow effect in the dissipative 
harmonic oscillator model, where results for the finite- 
temperature susceptibility turn out to be qualitatively 
incorrect in general. Applied to quantum phase transi- 
tions in bosonic impurity models, we have argued that 
the mass-flow effect introduces qualitative errors in the 
critical regime of mean- field quantum phase transitions, 
while it only leads to quantitative errors for interacting 
quantum criticality. A simple extension of the NRG algo- 
rithm allows to cure the mass-flow error asymptotically 
near the fixed points of interest. We have applied this 
modified algorithm to the sub-ohmic spin-boson model 
and found unambiguous signatures of mean-field behav- 
ior for s < 1/2, including a flow towards a Gaussian 
critical fixed point, Fig. 5b, and a susceptibility power 
law with mean-field exponent, Figs. 7b and 8b. We have 
thus resolved the discrepancy between results from NRG 
and those from other numerical methods. 10 ' 12 ~ 14 

As the conventional NRG is not capable of describing 
mean-field critical points, claims of non-mean-field be- 
havior in related 7 Ising-symmetric impurity models with 
sub-ohmic bosonic bath 5 ' 16 ' 17 need to be re-visited. 
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APPENDIX A: SPIN-BOSON MODEL: 
DETERMINING THE PROPER MASS-FLOW 
CORRECTION 

As discussed in Sec. V C, a general algorithmic solution 
to the mass-flow problem for interacting bosonic impu- 
rity models is not available. Instead, we have argued that 
an empirical correction via Eq. (12) within the NRG* al- 
gorithm is appropriate, with a prefactor n which depends 
on the fixed point of interest. 

Here we show the influence of k on the NRG* results for 
the sub-ohmic spin-boson model near criticality. Fig. 10 
displays NRG flows for a < a c (top) and local suscep- 
tibility data for various a ~ a c (bottom) for different 
values of k for s — 0.4. The central observation is that 
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FIG. 10: NRG* results for the spin-boson model with s = 0.4, A = 4, N b = 12, N a = 40, for different values of the flow- 
correction parameter k, Eq. (12). Top: Flow diagrams for a < ct c , specifically a — 0.3559, 0.3569, 0.364, 0.372 from a) to d). 
Bottom: Local susceptibilities % f° r various a ~ a c , evaluated with f3 — 1. The dashed (dash-dotted) lines show power laws 
with T -0,4 (T" 1/2 ) as reference. The a values are a) 0.355, 0.3559, 0.3559076, 0.356, b) 0.356, 0.3568, 0.356905, 0.3569096, 
0.35692, c) 0.36, 0.3635, 0.364, 0.365, d) 0.365, 0.369, 0.372, 0.374. Signatures of overcompensation, k > /to, are obvious in c) 
and d), while the critical power law \ ^ T~ 0A in a) implies undercompensation, k < kq. We conclude kq ~ 0.5. 



the behavior qualitatively changes when k is varied from 
0.3 to 1.0. 

In a) both the critical and delocalized NRG fixed 
points are clearly visible in the flow, and the critical 
x(T) oc T~ s . In b), the flow in the critical regime dis- 
plays a decreasing level spacing with increasing n, and 
no asymptotic T~ s power law is observed. Panels c) and 
d) show clear signs of overcompensation as discussed in 
Sec. VC, i.e., a non-monotonic flow (critical-localized- 
delocalized) and a corresponding non-monotonic x(T) 
near a c . Here, a precise determination of a c is impos- 
sible, and no critical power law in x(^) emerges. A 
detailed analysis of case a) shows that the behavior is 
qualitatively similar to that of the standard NRG. For 
a Gaussian fixed point, this would imply undercompen- 
sation. Together with the discussion in Sec. VC, these 
observations strongly suggest that the critical fixed point 
of the spin-boson for s = 0.4 is Gaussian, with kq w 0.5, 
see also the data in Fig. 5. Indeed, the critical x(T) m 
Fig. 10b does not follow x oc T~ x with x = 0.4 at the low- 
est T shown, but instead crosses over to larger x. A sim- 
ilar procedure for other s < 1/2 yields ko(s = 0.2) w 0.2 
and Ko(s = 0.3) » 0.35. 

In contrast, data for s = 0.6, Fig. 11, do not display 
a qualitative change when k is varied from 0.5 to 2.0. 
(k = data from the standard NRG are in Fig. 6.) 
Instead, for all K a stable critical fixed-point spectrum 
emerges, and the critical susceptibility follows x <x T~ s . 
This implies an interacting critical fixed point in the spin- 



a)K = 0.5 b)K = 2.0 
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FIG. 11: As in Fig. 10, but for s = 0.6. Here, no qualitative 
changes occur upon variation of k. The a values are 0.671007 
and 0.693268 in the top panel, and a) 0.67, 0.671, 0.67100924, 
0.67102, b) 0.69, 0.6932, 0.35739, 0.6932712, 0.6933 in the 
bottom panel. The dashed line shows a power law with J 1-0 - 6 . 



boson model for s — 0.6. An extrapolation of kq(s) sug- 



12 



gests /to(s = 0.6) i=a 0.8 (which, however, is not very accu- 
rate, see Sec. VC). We note that Figs. 11a and b differ 
quantitatively: The level structure at the critical fixed 
point is somewhat shifted, and the prefactor of the crit- 
ical power law of x(^) is 40% larger in b). This trend 
simply reflects that larger n reduces the order-parameter 
mass along the flow trajectory. 



APPENDIX B: FERMIONIC RESONANT-LEVEL 
MODEL 

The mass-flow effect is in principle also present in 
fermionic impurity models if the bath is particle-hole 
asymmetric in the low-energy limit. Consider the spinlcss 
resonant-level model 

Hrlm = e// t / + ^A 4 (/ t c 4 + h.c.) + ^w l c J t c l 

i i 

(Bl) 

with the bare level energy e/. As in Sec. II, one can de- 
fine a bath spectral density which, however, now 
generically has contributions at both positive and nega- 
tive frequencies. 

The solution for the / (impurity) Green's function is 



low-energy asymmetry may be quantified by looking at 



G / ( W ) = «/;/t)) 



i0+ - e f -T(lu) 



(B2) 



The impurity properties of this model are non-singular 
except at resonance, a = a c , where e/ + Rer(w = 0) = 0, 
i.e., where the renormalized / level coincides with the 
Fermi level. The properties near resonance have been 
studied extensively in Refs. 19,20 for the particle-hole 
symmetric case. 

A mass-flow error arises only for bath spectra J{u>) 
which arc particle-hole asymmetric at low energies. The 



a(iv) = 



J{lo) - J{-lo) 
J(u>) + J{-u>) 



(B3) 



A finite a(w^0) implies particle-hole asymmetry in lead- 
ing order. Otherwise the mass-flow error vanishes in the 
low-energy limit, this applies e.g. to a metallic fermionic 
bath spectrum with different positive and negative band 
cutoff energies. 

We have studied the mass-flow error for the resonant- 
level model (Bl) with a maximally particle-hole asym- 
metric power- law bath, i.e., J(u> > 0) oc lu s and J(u> < 
0) = 0. In analogy with Sec. IV, we then expect a large 
mass-flow error at resonance. Indeed, the resonance po- 
sition is shifted in a similar fashion as for the harmonic 
oscillator in Sec. IV. However, when comparing observ- 
ables, like the / level occupancy or its susceptibility, 
calculated at T n for the semi-infinite and the truncated 
chains with fixed ef, we find that the differences are tiny 
(less than 10~ 2 ), in stark contrast to the bosonic case. 
The reason for the small mass-flow error is rooted in both 
the character of the observables and the statistics of the 
particles. First, in the fermionic case all observables are 
related to two-particle propagators, in contrast to the 
bosonic \ 01 Eq. (8). Hence, the real part of T{u) never 
shows up as directly as in x, due to a convolution inte- 
gral. Second, the response of fermions at resonance is less 
singular than that of bosons. Therefore, deviations from 
the exact resonance condition have less consequences as 
compared to the bosonic case. 

In summary, the mass-flow error of NRG is present for 
particle-hole asymmetric fermionic problems as well, but 
practically has little effect for the observables we have 
checked. 
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